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ABSTRACT 


Risk/hazard  rating  methods  for  mountain  pine  beetle  (MPB)  in  ponderosa  pine  stands  were  compared, 
based  upon  number  of  stands  and  acres  that  rated  as  low,  medium  or  high  risk/hazard,  for  four  recent 
project  areas  on  the  Black  Hills  National  Forest.  Methods  compared  were  the  (1)  Stevens  model  described 
by  Stevens  et  al.  (1980)  except  that  the  basal  area  threshold  for  high  susceptibility  was  reduced  from  150 
to  120  based  upon  results  of  Schmid  and  Mata  (1992),  (2)  FORPLAN  model  used  in  the  Draft  Revised 
Black  Hills  National  Forest  Plan  (same  as  Stevens  model  except  for  assumption  of  only  one-stoned  stand 
structure),  (3)  PONBUG  model  for  the  Black  Hills  National  Forest  (same  as  Stevens  model  except  that 
data  for  trees  >  1  inch  in  diameter  at  breast  height  (DBH)  were  used  rather  >  5  inches  DBH),  (4)  Schmid 
model  described  by  Schmid  et  al.  (1994)  with  assumptions  added  to  extend  it  to  stands  with  quadratic 
mean  diameter  (QMD)  of  <  7  inches  and  >  13  inches,  and  (5)  Schmid  5  model  (same  as  Schmid  model 
but  recalculated  using  data  for  trees  >  5  inches  DBH  rather  than  >  1  inch  DBH). 

Considerable  differences  are  evident  in  results  of  the  five  methods.  The  FORPLAN  model  most  closely 
resembles  the  Stevens  model  although  some  stand  MPB  risk/hazard  ratings  shifted  from  low  to  medium 
and  from  medium  to  high.  In  general,  both  stands  and  acres  tended  to  increase  substantially  in  the  low 
MPB  risk/hazard  category  (frequently  doubling)  using  the  PONBUG  model  compared  to  the  Stevens 
model,  while  stands  and  acres  in  medium  and  high  MPB  risk/hazard  categories  tended  to  decline 
However,  the  pattern  was  somewhat  different  for  the  Tollgate  project  where  many  stands  had  an 
abundance  of  trees  <  5  inches  DBH.  The  Schmid  model  tended  to  shift  even  more  stands  and  acres  to  the 
low  MPB  risk/hazard  category  than  did  the  PONBUG  model.  In  general,  stands  and  acres  in  the  medium 
and  high  MPB  risk/hazard  categories  declined  for  the  Schmid  model  relative  to  the  Stevens  model  and  its 
variants,  with  a  few  exceptions,  most  notably  the  high  MPB  risk/hazard  category  in  the  Tollgate  project 
The  Schmid  5  model  tended  to  shift  more  stands  to  low  risk/hazard  or  high  risk/hazard  depending  upon 
project  area  considered  when  compared  to  the  other  models.  There  was  little  consistency  between  stands 
that  rated  high  in  the  Schmid  5  model  and  stands  that  rated  high  in  the  other  models. 

Differences  in  model  assumptions,  design,  and  data  parameters  affected  MPB  nsk/hazard  rating 
outcomes.  Assuming  only  one-storied  stands  (FORPLAN  model)  increased  risk/hazard  ratings  compared 
to  a  similar  (Stevens)  model  that  also  included  consideration  of  two-stoned  stands.  Model  results  differed 
when  using  data  for  trees  >  1  inch  DBH  versus  for  trees  >  5  inches  DBH  and  the  model  outcomes 
appeared  to  be  affected  by  abundance  of  trees  of  small  diameter  The  Schmid  model  in  particular  may  not 
adequately  reflect  MPB  risk/hazard  levels  for  uneven-aged  stands  and  those  with  many  trees  of  small 
DBH.  There  appear  to  be  inherent  differences  in  the  way  the  Stevens  model  and  its  variants  calculate 
MPB  risk/hazard  relative  to  the  Schmid  model  and  its  vanant  that  are  independent  of  which  data 
parameters  are  used  for  QMD  and  BA.  The  analyses  of  model  differences  presented  in  this  report  cannot 
identify  which  is  closest  to  correct  as  related  to  MPB  behavior  in  field  settings 
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INTRODUCTION 


Stevens  et  al.  (1980)  described  a  system  for  rating  ponderosa  pine  stands  in  the  Black  Hills  for 
susceptibility  to  mountain  pine  beetle  (MPB),  Dendroctonus  ponderosae  Hopkins,  using  parameters  of 
stand  structure  (one-  or  two-storied),  average  diameter  at  breast  height  (DBH),  and  average  basal  area 
(BA).  This  MPB  risk-rating  system  was  incorporated  into  a  model  subroutine  called  PONBUG,  which  is 
linked  to  the  Rocky  Mountain  Resource  Inventory  System  (RMRIS)  database  used  by  National  Forests  in 
the  Rocky  Mountain  Region,  so  that  MPB  risk  could  be  automatically  calculated  in  the  program  called 
RMSTAND.  Results  can  be  displayed  in  the  optional  field  “BUGRISK”  in  the  table  R2RIS  SITE 
Values  of  1,  3,  and  5  indicate  low,  medium,  and  high  MPB  stand  risk,  respectively.  The  program  was 
modified  in  1993  for  the  Black  Hills  and  Nebraska  National  Forests  based  upon  results  of  Schmid  and 
Mata  (1992).  They  suggested  that  120  sq.  ft.  per  acre  may  be  the  critical  BA  for  separating  high  risk 
stands  of  susceptible-sized  trees  in  the  Black  Hills,  rather  than  150  sq.  ft.  per  acre.  This  risk-rating 
program  is  often  used  by  district  personnel  during  project  planning  to  help  prioritize  pine  stands  for 
thinning  or  harvesting. 

As  Rapid  City  Service  Center  staff  became  involved  in  assisting  district  personnel  with  analyses  of  MPB 
risk,  it  became  apparent  that  the  PONBUG  model  did  not  always  produce  the  same  results  as  hand 
calculations  of  the  method  published  by  Stevens  et  al.  (1980).  Closer  examination  of  the  computer 
program  PONBUG  revealed  that  it  used  average  diameter  and  BA  of  trees  >  1  inch  QMD  (quadratic  mean 
diameter,  an  estimate  of  average  stand  diameter)  rather  than  those  for  trees  >  5  inches  in  diameter,  as 
specified  by  Stevens  et  al.  (1980).  In  addition,  Stevens  et  al.  (1980)  did  not  define  one-stoned  and  two- 
storied  stand  structure  in  their  rating  system.  In  order  to  use  the  stand  structure  parameter  in  the 
PONBUG  model,  a  stand  was  designated  as  two-storied  if  BA  of  trees  >  1  inch  QMD  was  >  30  sq.  ft.  per 
acre  and  number  of  live  trees  between  1  and  5  inches  QMD  was  >  250  per  acre.  Stands  without  these 
conditions  were  designated  as  one-stoned.  By  this  method,  most  stands  in  the  Black  Hills  would  be 
designated  as  one-storied.  Whether  or  not  this  is  an  accurate  portrayal  of  stand  structures  of  ponderosa 
pine  in  the  Black  Hills  is  debatable.  The  PONBUG  designations  for  stand  structure  may  not  represent 
what  Stevens  et  al.  (1980)  envisioned  when  they  described  their  rating  system.  The  omission  of  trees 
smaller  than  5  inches  DBH  from  their  model  suggests  that  they  considered  such  small  trees  to  be 
unimportant  for  risk  rating  mountain  pine  beetle  susceptibility  and  may  also  reflect  a  scarcity  of  such  trees 
in  stands  of  high  risk.  Their  two-storied  stands  were  likely  to  contain  scattered  mature  trees  with  an 
understory  of  pole-sized  trees,  typical  of  stands  found  some  years  following  shelterwood  or  seed  tree 
harvesting.  A  quantitative  means  of  describing  such  stands  for  use  in  risk/hazard  rating  systems  for 
mountain  pine  beetle  in  ponderosa  pine  has  not  been  developed  to  date. 

An  attempt  was  made  to  use  the  model  of  Stevens  et  al.  (1980)  as  modified  by  the  120  BA  change  based 
upon  results  of  Schmid  and  Mata  (1992)  in  analyses  of  alternatives  for  the  Proposed  Revised  Land  and 
Resource  Management  Plan  and  Draft  Environmental  Impact  Statement  (HIS)  for  the  Black  Hills  National 
Forest  (1994).  The  FORPLAN  model,  used  to  simulate  results  of  various  management  alternatives,  would 
not  handle  a  parameter  for  stand  structure  as  described  by  the  PONBUG  model,  so  an  assumption  was 
made  that  all  stands  were  one-storied  This  was  expected  to  shift  some  stands  to  a  higher  risk  rating  than 
would  be  produced  with  the  BA-modified  Stevens  model 

Several  months  after  the  Draft  Revised  Plan  and  EIS  were  released  for  public  comment,  Schmid  et  al. 
(1994)  published  an  alternative  method  of  rating  stands  for  MPB  hazard1  in  the  Black  Hills.  A  decision 


1  The  terms  ‘risk’  and  ‘hazard’  and  their  meanings  are  often  interchanged  and  confused  in  the  literature. 
As  used  in  this  report,  the  terms  both  refer  to  the  degree  to  which  stand  characteristics  provide  suitable 
habitat  for  development  of  MPB.  The  presence  or  absence  of  a  MPB  population  is  not  considered  for 
rating  determination.  Therefore,  risk/hazard  rating  as  used  herein  does  not  predict  when  or  how 
frequently  a  given  stand  will  be  attacked  by  MPB  Terminology  usage  of  rating  system  authors  was 
maintained  in  this  report;  i.e.,  Stevens  et  al.  (1980)  used  ‘risk’  whereas  Schmid  et  al.  (1994)  used 
‘hazard.’ 
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was  made  to  use  this  recent  research  in  analyses  for  the  Final  Rev  ised  Plan  and  EIS.  The  method  as 
described  by  Schmid  et  al.  (1994)  only  pertains  to  stands  of  7  to  13  inches  mean  DBH,  so  some 
assumptions  had  to  be  made  for  incorporation  in  the  SPECTRUM  model  (a  recent  replacement  for  the 
FORPLAN  model).  It  was  assumed  that  MPB  hazard  would  be  low  for  all  ponderosa  pine  stands  of  <  7 
inches  QMD.  An  assumption  also  was  made  that  stands  with  QMD  >  13  inches  would  have  similar 
hazard  ratings  as  for  stands  in  the  7  to  13  inch  DBH  range  (i.e.,  the  model  was  extended),  although  there 
is  some  evidence  to  suggest  that  MPB  nsk/hazard  may  decline  for  stands  with  larger  mean  diameters 
(Lessard  1982).  Data  are  limited  to  support  this  assertion  because  few  stands  with  QMD  >  13  inches  are 
present  in  the  Black  Hills.  Because  of  the  limited  acreage  of  large  diameter  ponderosa  pines,  the  overall 
hazard  ratings  for  forest  planning  purposes  was  not  expected  to  be  affected  much  even  if  the  assumption 
was  incorrect. 

The  release  of  the  Schmid  et  al.  (1994)  method  between  publication  of  the  draft  and  final  versions  of  the 
Revised  Forest  Plan  for  the  Black  Hills  National  Forest  raised  the  question  of  how  much  nsk/hazard 
ratings  from  the  Stevens  model  and  its  variations  (PONBUG,  FORPLAN)  might  differ  from  the 
nsk/hazard  ratings  of  the  Schmid  method.  It  also  seemed  appropriate  to  evaluate  the  methods  and  model 
assumptions  before  recommending  any  changes  in  the  PONBUG  program  used  for  project  analyses.  The 
model  of  Stevens  et  al.  (1980)  was  based  upon  relationships  reported  by  Sartwell  and  Stevens  (1975)  and 
general  observations  of  pine  stands  in  the  Black  Hills.  In  contrast,  Schmid  et  al.  (1994)  based  their  model 
on  sets  of  plots  of  even-aged  stands  comprised  of  susceptible-sized  trees  that  had  been  exposed  to  epidemic 
MPB  infestations.  This  report  compares  and  evaluates  the  risk/hazard  ratings  for  stands  in  four  project 
areas  as  derived  from  each  of  five  methods. 


METHODS 

Inventory  summary  data  generated  by  RMSTAND  using  the  RMRIS  database  were  obtained  for  the 
Soholt,  Thumb,  Kirk,  and  Tollgate  project  areas  (Fig.  1).  Data  on  stand  location  and  site  number,  number 
of  acres,  cover  type,  BA  and  QMD  of  trees  >  1  inch  in  diameter,  BA  and  QMD  of  trees  >  5  inches  in 
diameter,  and  numbers  of  live  trees  per  acre  for  trees  >  1  inch  and  >  5  inches  in  diameter  were  transferred 
from  the  Data  General  minicomputer  system  to  a  Microsoft  Excel  spreadsheet  on  a  desktop 
microcomputer.  Number  of  live  trees  per  acre  between  1  and  5  inches  in  diameter  was  calculated  by 
subtracting  the  number  of  trees  per  acre  >  5  inches  in  diameter  from  the  number  of  live  trees  per  acre  >  1 
inch  in  diameter,  and  then  stored  in  a  column  in  the  spreadsheet.  For  each  of  the  four  project  areas,  MPB 
risk/hazard  was  calculated  for  each  stand  of  ponderosa  pine  cover  type  that  had  complete  inventory 
information  on  BA,  QMD,  and  number  of  trees  per  acre  For  each  stand,  five  risk/hazard  rating  methods 
were  applied  and  compared  as  described  below. 

1.  Stevens  model — The  risk  rating  system  described  by  Stevens  et  al.  (1980)  was  used  with  the  exception 
that  the  high  susceptibility  threshold  for  BA  was  lowered  from  150  to  120  based  on  results  of  Schmid  and 
Mata  (1992).  Stand  structure  was  classified  as  two-storied  if  BA  of  trees  >  1  inch  DBH  was  >  30  sq.  ft 
per  acre  and  number  of  live  trees  between  1  and  5  inches  DBH  was  >  250  per  acre,  and  otherwise 
classified  as  one-storied  (same  as  programmed  in  PONBUG),  because  no  other  means  of  differentiating 
one-storied  and  two-storied  stands  of  ponderosa  pine  for  MPB  risk/hazard  rating  purposes  has  been 
described  or  used.  Data  for  mean  tree  diameter  and  basal  area  for  trees  >  5  inches  DBH  (designated  as 
QMD  5  and  BA_5,  respectively)  were  used  in  determination  of  risk  classes  for  mean  stand  diameter  and 
mean  stand  density.  Risk  classes  were  assigned  to  stand  characteristics  as  shown  in  Table  1.  Values  of 
risk  class  for  each  stand  characteristic  were  multiplied  together  to  obtain  a  stand  risk  value.  A  stand  risk 
rating  was  then  assigned  depending  upon  stand  risk  value  as  shown  in  Table  2.  Example:  A  one-stoned 
stand  (risk  class  =  3)  with  QMD  5  =  8  (risk  class  =  2)  and  a  BA_5  =  180  (risk  class  =  3)  would  have  a 
stand  risk  value  of  3  x  2  x  3  =  18.  This  would  just  make  the  high  stand  risk  rating  (code  5). 
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BLACK  HILLS 


Fig.  1.  Locations  of  four  projects  areas  (Soholt,  Thumb,  Kirk  and  Tollgate)  analyzed  for  mountain  pine 
beetle  risk/hazard  on  the  Black  Hills  National  Forest. 


Table  1.  Risk  classes  for  mountain  pine  beetle  nsk  rating  of 
ponderosa  pine  stands  in  the  Black  Hills  as  described  by  Stevens  et  al 
(1980)  and  modified  based  on  results  of  Schmid  and  Mata  (1992). 


1  =  low 

2  =  medium 

3  =  high 

Stand  structure 

— 

two-storied 

one-storied 

Mean  stand  DBH  (in.) 

<6 

6-10 

>10 

Stand  density  (BA) 

<80 

80-120 

>120 

Table  2.  Mountain  pine  beetle  stand  risk 
ratings  as  assigned  to  calculated  stand 
nsk  values  in  Black  Hills  ponderosa  pine 
(from  Stevens  et  al.  1980). 


Stand  risk  rating 

Stand  risk  value 

Low  =  1 

2-6 

Medium  =  3 

8-12 

High  =  5 

18-27 

2.  FORPLAN  model — Stand  structure  was  assumed  to  be  one-stoned  for  all  stands  (i.e.,  nsk  class  =  3). 
Otherwise,  calculations  were  the  same  as  for  the  Stevens  model  (including  use  of  QMD  5  and  BA_5  and 
the  values  in  Table  1). 

3.  PONBUG  model — MPB  risk/hazard  ratings  were  calculated  identically  to  the  Stevens  model  (using 
values  in  Table  1),  except  that  data  for  mean  tree  diameter  and  basal  area  for  trees  >  1  inch  DBH  were 
used  (designated  as  QMDl  and  BA_1)  instead  of  for  trees  >  5  inches  DBH. 

4.  Schmid  model — The  hazard  rating  method  described  by  Schmid  et  al.  (1994)  was  adapted  as  follows. 
Data  for  mean  tree  diameter  and  basal  area  for  trees  >  1  inch  DBH  was  used  (i.e.,  QMD  l  and  BA_1). 

For  QMD  l  of  <  7  inches,  stands  were  rated  as  low  hazard  (code  =1).  MPB  hazard  also  rated  low  (code 
=  1)  if  QMD  l  was  >  7  inches  and  growing  stock  level  (GSL)  was  <  80.  MPB  hazard  rated  medium 
(code  =  3)  if  QMD  l  was  >  7  inches  and  GSL  was  >  80  and  <  120.  MPB  hazard  was  rated  high  (code  = 
5)  if  QMD  l  was  >  7  inches  and  GSL  was  >  120.  For  stands  with  QMD  l  >  10  inches,  BA_1  of  live 
trees  was  directly  substituted  for  GSL.  For  stands  with  QMD  l  of  >  7  inches  and  <  10  inches,  GSL  was 
calculated  using  the  following  formula: 

GSL  =  (80  x  BA_  1)  -s-  [(7.76226  x  QMD_  1)  +  (0.85289  x  QMD_\2)  -  (0.07952  x  QMD _\ 3)  -  3.45624] 
(Carl  Edminster,  personal  communication). 

5.  Schmid  5  model — MPB  risk/hazard  ratings  were  calculated  identically  to  the  Schmid  model,  except 
that  data  for  mean  tree  diameter  and  basal  area  for  trees  >  5  inches  DBH  was  used  (i.e.,  QMD  5  and 
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BA_5)  rather  than  for  trees  >  1  inch  DBH.  This  method  was  intended  to  allow  more  direct  comparison  of 
results  for  the  Stevens  and  Schmid  methods  by  basing  calculations  on  the  same  data  parameters 

MPB  risk/hazard  ratings  for  the  PONBUG  and  Schmid  models  were  overlain  onto  ponderosa  pine 
stocking  charts  for  the  Black  Hills  to  illustrate  which  stand  conditions  fall  into  the  low,  medium,  and  high 
MPB  risk/hazard  categories.  Similar  charts  could  not  be  generated  for  the  Stevens,  FORPLAN,  and 
Schmid  5  models  because  stocking  charts  based  upon  data  for  trees  >  5  inches  DBH  are  not  available 

MPB  risk/hazard  rating  information  for  each  model  was  summarized  in  terms  of  numbers  and  percent  of 
stands  in  low,  medium,  and  high  MPB  risk/hazard  for  each  project  area  analyzed.  Additionally,  results 
for  each  model  were  summarized  by  number  and  percent  of  acres  rated  as  low,  medium,  and  high  in  each 
project  area. 

Results  of  MPB  risk/hazard  rating  calculations  were  transferred  to  ARC  VIE  W,  a  geographic  information 
system  (GIS)  program,  to  generate  maps  displaying  low,  medium,  and  high  risk  stands  using  the  various 
MPB  risk/hazard  rating  models  for  each  of  the  project  areas  analyzed 

RESULTS  AND  DISCUSSION 

MPB  risk  rating  categories  as  produced  by  PONBUG  are  displayed  on  stocking  charts  for  Black  Hills 
ponderosa  pine  for  one-stoned  and  two-storied  stands  (Figs.  2  and  3).  Hazard  rating  categones  of 
Schmid  et  al.  (1994)  including  added  assumptions  for  forest  planning  purposes  are  illustrated  on  the 
stocking  chart  in  Fig.  4.  The  Schmid  model  has  a  narrower  range  of  stand  conditions  that  fall  into  the 
medium  MPB  nsk/hazard  category.  The  range  of  high  MPB  risk/hazard  stand  conditions  is  somewhat 
smaller  for  the  Schmid  model  than  for  one-storied  stands  using  PONBUG,  but  is  considerably  larger  than 
that  described  for  two-storied  stands  using  PONBUG.  The  range  of  stand  conditions  rated  as  low  by  the 
Schmid  model  is  larger  than  the  ranges  encompassed  by  low  ratings  for  both  one-storied  and  two-stoned 
stands  in  PONBUG 

Considerable  differences  in  the  number  and  percent  of  stands  (Table  3)  and  number  and  percent  of  acres 
(Table  4)  rating  as  low,  medium,  and  high  risk/hazard  for  mountain  pine  beetle  are  evident  for  the  five 
methods  evaluated  in  each  of  the  four  project  areas.  Results  of  the  FORPLAN  model  most  closely 
resemble  those  of  the  Stevens  model,  with  some  stand  ratings  shifting  from  low  to  medium  and  from 
medium  to  high  (Tables  3-5).  None  of  the  individual  stand  ratings  shifted  to  lower  values  when 
comparing  results  for  the  FORPLAN  model  to  the  Stevens  model.  This  general  shift  to  higher  MPB 
nsk/hazard  ratings  was  expected  based  upon  the  assumption  of  only  single-stoned  stands  for  the 
FORPLAN  model.  The  FORPLAN  model,  apparently  because  of  this  assumption,  generally  produced  the 
greatest  number  of  stands  that  rated  high  risk/hazard  of  the  five  models,  except  in  the  Tollgate  project 
where  it  had  the  second  highest  number  of  stands  rated  as  high  risk/hazard  (Table  3).  Departures  of  the 
FORPLAN  model  from  the  Stevens  model  were  greatest  for  the  Tollgate  and  Kirk  projects,  located  in  the 
northern  Black  Hills  where  growing  conditions  are  better  for  ponderosa  pine. 

Both  stands  and  acres  tended  to  increase  substantially  in  the  low  MPB  risk/hazard  category  (frequently 
doubling)  using  the  PONBUG  model  compared  to  the  Stevens  model,  on  which  the  PONBUG  model 
supposedly  was  based  (Tables  3-4).  Simultaneously,  stands  and  acres  in  medium  and  high  MPB 
risk/hazard  categories  tended  to  decline.  These  relations  held  for  the  Soholt,  Thumb,  and  Kirk  projects; 
however,  stands  and  acres  rated  as  low  and  high  MPB  risk/hazard  declined  and  medium  MPB 
nsk/hazard  increased  in  the  Tollgate  project  area  using  the  PONBUG  model  compared  to  the  Stevens 
model  (Tables  3-4).  Changes  in  individual  stand  ratings  both  increased  and  decreased  for  the  PONBUG 
model  relative  to  the  Stevens  model;  however,  a  majority  of  stands  that  changed  ratings  went  from 
medium  to  low  MPB  risk/hazard  in  the  Soholt,  Thumb,  and  Kirk  projects  while  the  greatest  change  was 
from  low  to  medium  MPB  risk/hazard  in  the  Tollgate  project  (Table  5).  Differences  between  the  Stevens 
and  PONBUG  models  would  be  the  result  of  using  inventory  data  based  on  trees  >  1  inch  DBH  in  the 
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Black  Hills  displayed  on  a  timber  stocking  chart. 


8 


Fig.  3.  Mountain  pine  beetle  risk  categories  as  produced  by  the  PONBUG  model  for  two-storied  stands  of  ponderosa  pine  in  the  Black 
Hills  displayed  on  a  timber  stocking  chart 


9 


with  added  assumptions  and  displayed  on  a  timber  stocking  chart. 


Table  3.  Number  and  percent  of  ponderosa  pine  stands  in  each  MPB  nsk/hazard  rating 
category  as  determined  by  five  calculation  methods  in  four  project  areas 


PROJECT  AREA 

NUMBER  (%)* 

OF  STANDS 

LOW 

MED 

HIGH 

TOTAL 

SOHOLT 

Stevens  Model 

6  (22) 

12  (44) 

9  (33) 

27 

FORPLAN  Model 

6  (22) 

11  (41) 

10  (37) 

27 

PONBUG  Model 

10  (37) 

12  (44) 

5  (19) 

27 

Schmid  Model 

13  (48) 

12  (44) 

2  (  7) 

27 

Schmid  5  Model 

15  (56) 

9  (33) 

3  (11) 

27 

THUMB 

Stevens  Model 

78  (31) 

154  (60) 

23  (  9) 

255 

FORPLAN  Model 

76  (30) 

156  (61) 

23  (  9) 

255 

PONBUG  Model 

153  (60) 

90  (35) 

12  (  5) 

255 

Schmid  Model 

196  (77) 

55  (22) 

4  (  2) 

255 

Schmid  5  Model 

208  (82) 

44  (17) 

3  (  1) 

255 

KIRK 

Stevens  Model 

39  (15) 

149  (56) 

78  (29) 

266 

FORPLAN  Model 

33  (12) 

128  (48) 

105  (39) 

266 

PONBUG  Model 

79  (30) 

119  (45) 

68  (25) 

266 

Schmid  Model 

148  (56) 

55  (21) 

63  (24) 

266 

Schmid  5  Model 

72  (27) 

114  (43) 

80  (30) 

266 

TOLLGATE 

Stevens  Model 

88  (35) 

132  (53) 

28  (11) 

248 

FORPLAN  Model 

45  (18) 

111  (45) 

92  (37) 

248 

PONBUG  Model 

73  (29) 

153  (62) 

22  (  9) 

248 

Schmid  Model 

95  (38) 

82  (33) 

71  (29) 

248 

Schmid  5  Model 

58  (23) 

78  (33) 

112  (45) 

248 

*  Percentages  of  low,  medium,  and  high  may  not  total  100  due  to  rounding. 
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Table  4.  Number  and  percent  of  acres  of  ponderosa  pine  in  each  MPB  risk/hazard  rating 
category  as  determined  by  five  calculation  methods  in  four  project  areas 


PROJECT  AREA 

NUMBER  (%) 

*  OF  ACRES 

LOW 

MED 

HIGH 

TOTAL 

SOHOLT 

Stevens  Model 

331  (15) 

1127  (51) 

757  (34) 

2215 

FORPLAN  Model 

331  (15) 

1078  (49) 

806  (36) 

2215 

PONBUG  Model 

928  (42) 

988  (45) 

299  (13) 

2215 

Schmid  Model 

1085  (49) 

1037  (47) 

93  (  4) 

2215 

Schmid  5  Model 

1382  (62) 

706  (32) 

127  (6) 

2215 

THUMB 

Stevens  Model 

4008  (29) 

8565  (62) 

1305  (  9) 

13878 

FORPLAN  Model 

3918  (28) 

8655  (62) 

1305  (  9) 

13878 

PONBUG  Model 

8393  (60) 

4806  (35) 

679  (  5) 

13878 

Schmid  Model 

10518  (76) 

3209  (23) 

151  (  1) 

13878 

Schmid  5  Model 

11082  (80) 

2590  (19) 

206  (  1) 

13878 

KIRK 

Stevens  Model 

1625  (16) 

5813  (59) 

2481  (25) 

9919 

FORPLAN  Model 

1402  (14) 

5260  (53) 

3257  (33) 

9919 

PONBUG  Model 

3206  (32) 

4723  (48) 

1990  (20) 

9919 

Schmid  Model 

5627  (57) 

2106  (21) 

2186  (22) 

9919 

Schmid  5  Model 

2900  (29) 

4567  (46) 

2452  (25) 

9919 

TOLLGATE 

Stevens  Model 

4063  (45) 

4088  (45) 

978  (11) 

9129 

FORPLAN  Model 

1894  (21) 

4060  (44) 

3175  (35) 

9129 

PONBUG  Model 

3343  (37) 

5096  (56) 

690  (  8) 

9129 

Schmid  Model 

3987  (44) 

2899  (32) 

2243  (25) 

9129 

Schmid  5  Model 

2273  (25) 

3105  (34) 

3751  (41) 

9129 

*  Percentages  of  low,  medium,  and  high  may  not  total  100  due  to  rounding. 
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Table  5.  Number  of  ponderosa  pine  stands  changing  rating  category  for  mountain  pine  beetle  nsk/hazard 
between  specified  model  pairs  in  four  project  areas 


Project 

Area 

Model  Comparison 

Number  of  Stands  Changing  Risk/hazard  Rating  Category 

Increase  in  Risk/hazard 

Decrease  in  Risk/hazard 

Low  to 
Medium 

Low 

to 

High 

Medium 
to  High 

Medium 
to  Low 

High  to 
Medium 

High 
to  Low 

No  Change 

Soholt 

Stevens  to  FORPLAN 

0 

0 

1 

0 

0 

0 

26 

Stevens  to  PONBUG 

1 

0 

1 

5 

5 

0 

15 

Stevens  to  Schmid 

3 

0 

1 

9 

7 

1 

6 

PONBUG  to  Schmid 

4 

0 

0 

6 

2 

1 

14 

Schmid  to  Schmid  5 

2 

1 

0 

5 

0 

0 

19 

Stevens  to  Schmid  5 

1 

0 

2 

8 

6 

2 

8 

Thumb 

Stevens  to  FORPLAN 

2 

0 

0 

0 

0 

0 

253 

Stevens  to  PONBUG 

4 

0 

0 

79 

11 

0 

161 

Stevens  to  Schmid 

14 

0 

1 

131 

19 

1 

89 

PONBUG  to  Schmid 

16 

0 

3 

59 

11 

0 

166 

Schmid  to  Schmid  5 

2 

2 

0 

16 

3 

0 

232 

Stevens  to  Schmid  5 

8 

0 

2 

135 

19 

3 

88 

Kirk 

Stevens  to  FORPLAN 

6 

0 

26 

0 

0 

0 

234 

Stevens  to  PONBUG 

2 

0 

5 

43 

15 

0 

201 

Stevens  to  Schmid 

3 

0 

16 

103 

22 

9 

113 

PONBUG  to  Schmid 

4 

0 

19 

64 

16 

8 

155 

Schmid  to  Schmid  5 

48 

33 

0 

5 

17 

0 

163 

Stevens  to  Schmid  5 

0 

0 

35 

29 

30 

4 

168 

12 


Table  5  continued.  Number  of  ponderosa  pine  stands  changing  rating  category  for  mountain  pine  beetle 
nsk/hazard  between  specified  model  pairs  for  four  project  areas. 


Project 

Area 

Model  Comparison 

Number  of  Stands  Changing  Risk/hazard  Rating  Category 

Increase  in  Risk/hazard 

Decrease  in  Risk/hazard 

Low  to 
Medium 

Low 

to 

High 

Medium 
to  High 

Medium 
to  Low 

High  to 
Medium 

High 
to  Low 

No  Change 

Tollgate 

Stevens  to  FORPLAN 

43 

0 

64 

0 

0 

0 

141 

Stevens  to  PONBUG 

21 

0 

6 

5 

11 

1 

204 

Stevens  to  Schmid 

21 

0 

48 

27 

4 

1 

147 

PONBUG  to  Schmid 

5 

0 

52 

27 

3 

0 

161 

Schmid  to  Schmid  5 

29 

19 

31 

11 

9 

0 

149 

Stevens  to  Schmid  5 

22 

13 

75 

5 

4 

0 

129 

All 

Stevens  to  FORPLAN 

51 

0 

91 

0 

0 

0 

654 

Stevens  to  PONBUG 

28 

0 

12 

132 

42 

1 

581 

Stevens  to  Schmid 

41 

0 

66 

270 

52 

12 

355 

PONBUG  to  Schmid 

29 

0 

74 

156 

32 

9 

496 

Schmid  to  Schmid  5 

81 

55 

31 

37 

29 

0 

563 

Stevens  to  Schmid  5 

31 

13 

114 

177 

59 

9 

393 

13 


PONBUG  model  rather  than  using  comparable  data  based  upon  trees  >  5  inches  DBH.  Including  trees  of 
smaller  (1  -4.9  inch)  DBH  in  the  data  would  tend  to  reduce  mean  diameter  or  QMD  and  increase  BA 
values  used  in  MPB  nsk/hazard  rating  calculations  The  pattern  of  change  in  stand  ratings  likely  was 
different  for  the  Tollgate  project  compared  to  the  other  locations  due  to  very  high  numbers  of  small  trees 
in  many  stands.  Trees  per  acre  in  the  1  -  4.9  inch  diameter  class  averaged  1,781  per  stand  for  the  Tollgate 
project  compared  to  92,  58,  and  219  for  the  Soholt,  Thumb,  and  Kirk  projects,  respectively.  For  the 
Soholt,  Thumb,  and  Kirk  projects,  the  differences  in  MPB  risk/hazard  rating  for  these  two  methods  likely 
were  due  primarily  to  lower  QMDl  values  relative  to  QMD5  values  For  the  Tollgate  project,  there 
likely  was  a  more  complex  mix  of  lower  QMD  l  values  and  higher  BA_1  values  that  affected  results  of 
the  PONBUG  model  compared  to  the  QMD5  and  BA_5  values  used  by  the  Stevens  model. 

The  Schmid  model  tended  to  shift  even  more  stands  and  acres  to  the  low  MPB  risk/hazard  category  than 
the  PONBUG  model  did  relative  to  the  Stevens  model  (Tables  3-4).  In  general,  stands  and  acres  in  the 
medium  and  high  MPB  risk/hazard  categories  declined  for  the  Schmid  model  relative  to  the  Stevens, 
FORPLAN,  and  PONBUG  models,  with  a  few  exceptions,  most  notably  the  high  MPB  nsk/hazard 
category  in  the  Tollgate  project  (Tables  3-4).  Most  stands  that  changed  ratings  in  the  Schmid  model 
compared  to  the  Stevens  and  PONBUG  models  went  from  the  medium  to  low  nsk/hazard  categories  in  the 
Soholt,  Thumb  and  Kirk  projects,  but  from  medium  to  high  MPB  risk/hazard  in  the  Tollgate  project 
(Table  5).  The  Schmid  model  shifted  MPB  risk/hazard  ratings  relative  to  the  PONBUG  model  as  would 
generally  be  expected  from  examination  of  the  stocking  charts  in  Figs.  2-4,  given  a  representative  array  of 
stand  conditions.  This  suggests  that  the  stands  in  the  four  project  areas  analyzed  in  this  report  provided  a 
reasonable  variety  of  stand  conditions.  The  variation  in  direction  of  change  (increase  or  decrease)  for  the 
Schmid  model  relative  to  the  PONBUG  and  Stevens  models  for  the  high  MPB  risk/hazard  category  for 
different  project  areas  likely  reflects  the  absence  of  stand  structure  as  a  model  component  in  the  Schmid 
model.  The  high  numbers  of  small  (1  -  4.9  inch  DBH)  trees  in  the  Tollgate  project  resulted  in  86%  of  the 
stands  being  classified  as  two-stoned  in  the  Stevens  and  PONBUG  models  compared  to  7,  2,  and  18 
percent,  respectively,  for  the  Soholt,  Thumb,  and  Kirk  projects.  A  preponderance  of  one-storied  stands  (as 
found  in  the  Soholt,  Thumb,  and  Kirk  projects)  might  tend  to  produce  fewer  stands  with  high  MPB 
nsk/hazard  ratings  as  compared  to  a  method  such  as  the  Schmid  model  that  considers  only  average  tree 
diameters  and  densities.  Although  Schmid  et  al.  (1994)  did  not  specify-  a  minimum  tree  diameter  to  be 
used  in  model  calculations,  their  method  was  based  upon  data  from  stands  that  contained  few  trees  <  4 
inches  DBH.  The  method  is  geared  toward  even-aged  stands  where  most  trees  are  of  susceptible  size 
Therefore,  the  Schmid  model  may  not  adequately  reflect  MPB  risk/hazard  levels  for  uneven-aged  stands 
and  those  with  many  trees  of  small  DBH  The  assumption  added  to  classify  stands  with  QMD  <  7  inches 
as  being  low  risk/hazard  may  also  lead  to  more  low  ratings  than  might  be  reasonably  expected  For 
example,  a  stand  of  high  basal  area  with  trees  in  the  4  -  10  inch  DBH  classes  might  have  a  QMD  just 
below  7  inches  and  yet  have  many  suitable  host  trees  for  MPB  infestation.  Interestingly,  the  Schmid 
model  tends  to  provide  more  conservative  estimates  of  MPB  risk/hazard  in  ponderosa  pine  than  the 
Stevens  model  and  its  variations. 

No  consistent  pattern  was  found  in  how  risk/hazard  ratings  changed  if  small  trees  (<  5  inches)  were 
ignored,  as  in  the  Schmid  5  model  (Tables  3-5).  The  Schmid  5  model  tended  to  shift  more  stands  and 
acres  into  the  low  MPB  risk/hazard  category  for  the  Soholt  and  Thumb  project  areas,  but  more  into  the 
high  MPB  risk/hazard  category  for  the  Tollgate  project  relative  to  the  other  four  methods  (Tables  3-4). 

For  the  Kirk  project,  ratings  produced  by  the  Schmid  5  model  tended  to  be  intermediate  compared  to  the 
other  models  (Tables  3-4).  For  the  Soholt  and  Thumb  projects,  the  greatest  change  was  from  medium  to 
low  MPB  risk/hazard  when  comparing  both  the  Stevens  and  Schmid  models  to  the  Schmid  5  model  (Table 
5).  For  the  Tollgate  project,  substantially  more  stands  increased  in  MPB  risk/hazard  rating  than 
decreased  for  the  same  model  comparisons,  with  the  greatest  change  being  from  medium  to  high  MPB 
risk/hazard  (Table  5).  In  the  Kirk  project,  the  greatest  change  from  the  Schmid  model  to  the  Schmid  5 
model  was  from  low  to  medium  MPB  risk/hazard,  while  many  stands  also  went  from  low  to  high  MPB 
risk/hazard  (Table  5).  In  contrast,  the  greatest  change  from  the  Stevens  model  to  the  Schmid  5  model  in 
the  Kirk  project  was  from  medium  to  high  MPB  risk/hazard;  however,  more  stands  decreased  in  MPB 
nsk/hazard  than  increased  overall  (Table  5).  There  appears  to  be  little  relationship  between  which  stands 
rated  high  for  the  Schmid  5  model  relative  to  the  other  methods  (Figs.  5-24).  The  lack  of  correspondence 
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between  stand  risk/hazard  ratings  for  the  Schmid  5  model  relative  to  the  Stevens  model  suggests  that  there 
are  inherent  differences  in  the  way  the  Stevens  and  Schmid  models  calculate  MPB  nsk/hazard  that  are 
independent  of  which  data  parameters  for  QMD  and  BA  were  used. 

Unfortunately,  the  analyses  of  model  differences  presented  in  this  report  cannot  differentiate  which  is 
closest  to  correct  as  related  to  MPB  behavior  in  field  settings  Such  a  determination  requires  waiting  for 
natural  outbreaks  of  MPB  populations  in  Black  Hills  ponderosa  pine  stands  where  damage  levels  can  be 
compared  to  stand  risk/hazard  levels  predicted  by  the  v  arious  models. 

The  Schmid  model  was  used  for  analysis  of  effects  of  management  alternatives  in  the  Final  EIS  for  the 
1997  Revision  of  the  Land  and  Resource  Management  Plan  of  the  Black  Hill  National  Forest,  rather  than 
the  FORPLAN  model  that  was  used  in  the  Draft  EIS.  This  decision  was  made  based  upon  the  currency  of 
the  research  and  because  the  method  was  developed  from  data  obtained  from  stands  exposed  to  MPB 
infestation  rather  than  being  derived  mostly  from  general  observations.  As  long  as  the  same  method  was 
applied  to  all  alternatives  considered  in  the  plan,  the  relative  ranking  of  alternatives  should  be  valid, 
within  the  limits  of  our  knowledge,  even  if  the  actual  magnitude  of  acres  ranked  as  low,  medium,  and 
high  are  not  entirely  accurate.  The  PONBUG  model  has  continued  to  be  used  by  district  staffs  for 
conducting  environmental  analyses  of  proposed  project  areas  on  the  Black  Hills  National  Forest. 
Consideration  should  be  given  to  correcting  the  PONBUG  model  to  match  criteria  developed  by  Stevens  et 
al.  (1980)  or  to  switching  from  the  use  of  the  PONBUG  model  to  the  Schmid  model  for  project-level 
planning,  depending  on  which  of  the  models  appears  to  better  reflect  MPB  infestations  in  the  Black  Hills 
in  the  future.  Further  refinement  of  either  of  these  models  is  also  needed,  especially  as  management  of  the 
ponderosa  pine  type  shifts  from  application  of  predominately  even-aged  silvicultural  systems  with  even 
tree  spacing  to  greater  emphasis  on  variable  density  thinning  and  possibly  uneven-aged  harvesting  to 
enhance  structural  diversity  and  visual  objectives. 
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Fig.  5.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (G1S)  maps 
as  generated  by  the  Stevens  model  for  the  Soholt  project  area. 
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Fig.  6.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS)  maps 
as  generated  by  the  FORPLAN  model  for  the  Soholt  project  area. 


17 


B) 


MPB  Risk/hazard 
|  |  Unrated 

Low 
Medium 
High 


0.5 


0 


0.5 


1.5 


2.5  Kilometers 


Fig.  7.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS)  maps 
as  generated  by  the  PONBUG  model  for  the  Soholt  project  area. 
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Fig  8  Mountain  pine  beetle  nsk/hazard  ratings  displayed  on  geographic  information  system  (GIS)  maps 
as  generated  by  the  Schmid  model  for  the  Soholt  project  area. 
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Fig.  9.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS)  maps 
as  generated  by  the  Schmid  5  model  for  the  Soholt  project  area. 
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Fig.  10.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  Stevens  model  for  the  Thumb  project  area. 
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Fig.  1 1 .  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  FORPLAN  model  for  the  Thumb  project  area. 
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Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
generated  by  the  PONBUG  model  for  the  Thumb  project  area. 
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Fig  13.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  Schmid  model  for  the  Thumb  project  area. 
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Fig.  14  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
■  maps  as  generated  by  the  Schmid  5  model  for  the  Thumb  project  area 
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Fig  15.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  Stevens  model  for  the  Kirk  project  area. 
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Fig.  16.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  FORPLAN  model  for  the  Kirk  project  area 
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Fig.  17.  Mountain  pine  beetle  nsk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  PONBUG  model  for  the  Kirk  project  area. 
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Fig  18.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  Schmid  model  for  the  Kirk  project  area 


29 


Fig  19.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  Schmid  5  model  for  the  Kirk  project  area. 
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Fig.  20  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  Stevens  model  for  the  Tollgate  project  area 
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Fig.  21.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  FORPLAN  model  for  the  Tollgate  project  area 
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Fig  22.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  PONBUG  model  for  the  Tollgate  project  area. 
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Fig  23.  Mountain  pine  beetle  nsk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  Schmid  model  for  the  Tollgate  project  area. 
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Fig  24.  Mountain  pine  beetle  risk/hazard  ratings  displayed  on  geographic  information  system  (GIS) 
maps  as  generated  by  the  Schmid  5  model  for  the  Tollgate  project  area 
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